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"^ ' Abstract 

o : 

^^ ■ Results on recent QCD measurements performed at the Tevatron pp Collider 

at ^/s = 1.96 TeV are here reported. The inclusive jet and dijet mass cross 
S ■ sections are compared to NLO pQCD calculations and to Run I results. The 

production rates and kinematic properties of W + jets production processes are 
compared to "enhanced" LO theoretical predictions. Non-perturbative "soft" 
interactions leading to the underlying event are studied and compared to QCD 
Monte Carlo phenomenological models. 
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1 Introduction 

Measurements aimed to test the predictions of Quantum Chromodynamics 
(QCD), the currently accepted theory of the strong interactions among quarks 
and gluons, represent a very important part of the physics program carried out 
at the Tevatron pp Collider. The large amount of data expected to be accumu- 
lated during Run II and the increase in center-of-mass energy (-^/s) from 1.8 

to 1.96 TeV, give CDFLU' and D0 ' experiments an unique opportunity to 
make precision tests of next-to-leading order perturbative QCD (NLO pQCD) 
and, by looking for deviations from theory, to search for new particles and new 
interactions down to a distance scale of ~ 10~^^ m. 

An optimal understanding of QCD in hadron collisions allows to improve 
the constraints on the fundamental parameters of the theory, as and the parton 
distribution functions (PDFs) ; results in a better control on the standard QCD 
production which represents the main background for most of the processes 
of interest, such as top and Higgs production; gives phenomenological input 
for the modeling of the non-perturbative regime (where the theory fails in its 
predictivity), such as the "soft" interactions generating the underlying event 
which accompanies the "hard" collision. 

2 Inclusive Jet Cross Section 

2.1 Experience from Run I 

During Run I, the CDF and D0 Collaborations performed several QCD mea- 
surements which, in general, were found to be in reasonable agreement with 
theoretical expectations. However, initial inclusive jet cross section measure- 
ments showed an excess of data at high Et over NLO pQCD predictions, 
which raised great interest among the high energy physics community and 
stimulated a reevaluation of the uncertainties associated to theoretical calcu- 
lations ™ . Subsequent studies have demonstrated that such excess can be 
explained within the Standard Model in terms of a larger than expected gluon 
distribution at high x. A better agreement of data versus theory was actually 
observed in subsequent measurements involving an increased data sample when 
PDFs with an enhanced gluon contribution at high x (CTEQ4HJ) were consid- 
ered!^ '^ (see figm left). Given the large uncertainty in the NLO calculations 
arising from the flexibility allowed by current knowledge of PDFs as well as the 
uncertainty in the experimental results, CDF and D0 data were found to be 
consistent between them, with previous measurements and with NLO pQCD. 
Most recent global PDFs fits have used results from these Run I analyses 
so to include the Tevatron high Et jet data in the determination of the high 
X gluon distribution. In particular, by involving jets in a range of rapidity 
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Figure 1: Left: inclusive jet cross section as measured in CDF using 87 pb^^ 
of Run I data at ^/s = 1.8 TeVuj. Right: uncertainty hand for the gluon 
distribution function at Q^ = 10 GeV in the CTEQ6 M se t*^ ; comparisons 
to CTEQ5M1 (solid), CTEQ5HJ (dashed) and MRSTOl^ (dotted) are also 
shown in terms of ratios to CTEQ6M. 



intervals up to the forward region, the D0 measurement allowed to constrain 
the partons over a much wider x kinematical range. These new PDF sets 

(CTEQgI^ and MRST02I3) represent the most complete information available 
for Run II QCD predictions, but are still characterized by a big uncertainty on 
the gluon distribution at high x (see fign] right). 

With larger data samples collected with a higher cross section^ , as well as 
with the improved performances of the upgraded CDF and D0 detectors. Run 
II jet measurements are expected to reduce this uncertainty so to give the best 
constraint on PDFs before the LHC era. 



2.2 Preliminary Run II measurements 

Both CDF and D0 experiments have recently performed preliminary Run II 
measurements of the inclusive jet cross section. 

In CDF, jets were reconstructed using the same Run I cone algorithm 

(JetClu-', R = 0.7) in the rapidity range 0.1 < \r]\ < 0.7; the same Run I 



^The increase in ^/s from 1.8 to 1.96 TeV leads to a significant increase in 
jet cross section at high Pt"^ , about a factor of 5 at Pt'^ ~ 600 GeV. 
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Figure 2: Inclusive jet cross section as measured in CDF analyzing 177 pb^^ 
of Run II data at ^/s = 1.96 TeV. Upper left: the results are compared on 
a log scale to the NLO pQCD prediction from EKS (fi = Ej^/2, Rsep = 1-3, 
CTEQ6.1M PDFs). Upper right: ratio of data to theoretical prediction on 
a linear scale (the upper plot is a zoom obtained by excluding the last bin). 
Lower: Ratio of Run II inclusive jet cross section to Run I one compared to 
theoretical prediction. The increase in cross section due to the increase in ^/s 
is evident. Data points include the statistical error, the error band represents 
the change in the cross section due to the 3% energy scale shift corresponding 
to the dominant experimental uncertainty. The two solid lines represent the 
dominant theoretical uncertainty due to PDFs. 



correction procedures were also used to properly correct thejet spectrum for 
both detector and physics effects affecting jet measurements'^. Four different 
samples (from 20, 50, 70 and 100 GeV jet Et trigger threshold requirements). 



were combined after correcting for trigger prescale and efficiency effects giving 
a jet Et spectrum already extending the Run I reach by about 150 GeV and 
spanning approximately 9 orders of magnitude ^ . The results reported here and 
summarized in fig|21correspond to an integrated luminosity of 177 pb~^ collected 
during the period 2002-3. In general, a reasonably good agreement is found 
within errors between data and NLO pQCD theory (EKSl^) using the latest 

CTEQ PDF sets '^ . Systematic uncertainties are found to be dominantcd by 
the error on jet energy scale in data and on high x gluon distribution function 
in theoretical calculations. Work is in progress in order to reduce the energy 

scale uncertainty, to use other jet algorithms like MidPoint and Kt'"'^ and to 
extend the analysis to forward jets. 

A completely new analysis was performed in D0. Jets were reconstructed 
with an "optimized" cone algorithm (MidPoint '-'■i^ , R = 0.7) and then cor- 
rected with a new set of jet corrections as derived from 7-jet and dijet balancing 

studies on data'--" as well as from Monte Carlo simulations. The jet energy 
scale derived with this method was characterized by large statistical and sys- 
tematic uncertainties increasing with jet energy due to extrapolation (mostly 
generated by small 7-l-jet sample statistics above 200 GeV). Jet measurements 
were restricted to the central region (|?7| < 0.5) to limit the impact of these 
uncertainties which are the dominant ones, but are expected to be gradually 
reduced with an increasing amount of collected data. Finally, jet resolution 
effects were taken into account by unsmearing the measured jet spectrum with 
an unfolding procedure. Fig|21 summarizes the results of a preliminary mea- 
surement from the analysis of four different jets samples (from 25, 45, 65 and 
95 GeV jet Pt trigger thresholds) based on an integrated luminosity of 34 pb^^. 
An agreement within the rather large uncertainties (dominated by the system- 
atic error on the jet energy scale) is observed between data and NLO theory 

(JETRAdIEI)) using CTEQ6 and MRSTOl PDF sets. 

3 Dijet Invariant Mass Distribution 

The dijet mass cross section measurement represents another powerful test of 
NLO pQCD. It is also sensitive to new physics such as quark compositeness 
(expected to give departure from theory at high masses) and new models be- 
yond the Standard Model predicting new particles decaying to dijet (expected 
to generate a peak in the dijet invariant mass spectrum at high values). 

The CDF and D0 Collaborations have both performed a preliminary Run 
II measurement of the dijet mass spectrum using similar criteria as well as 
the same jet trigger samples as the ones involved for the inclusive jet cross 
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Et = 666 GeV is the highest jet Et so far observed. 
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Figure 3: Inclusive jet cross section as measured in _D0 analyzing 34 pb ^ of 
Run II data at ^/s — 1.96 TeV. Upper: the results are compared on a log 
scale to the NLO pQCD prediction from JETRAD (^jl = P™°^/2, Rsep = 1-3, 
CTEQ6M PDFs). Lower: linear comparison of data to theoretical prediction 
with CTEQ6M (left) and MRSTOl (right) PDF sets. Data points include the 
statistical error, the error band represents the total systematic uncertainty. 



section study. In particular, considering the two leading jets per event in order 
to reconstruct the invariant mass observable, the same Run I cone algorithm 
and analysis strategies (Jetclu R — 0.7 jets reconstructed up to \ri\ < 2.0 and 
corrected a the parton level) were used in CDF, while a completely new analysis 
was performed in D0 using the MidPoint algorithm and applying the same jet 
scale corrections and unsmearing procedures as the ones used for the inclusive 
jet cross section measurement to jets reconstructed in the same rapidity range 
i\v\ < 0.5). 

Results from CDF, corresponding to an integrated luminosity of 75 pb ^, 
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Figure 4: Upper: comparison of the dijet mass spectrum as measured in CDF 
using 106 pb^^ and 75 pb^^ of Run I and Run II data respectively (left): the 
CDF Run II over Run I dijet mass cross section ratio is compared to a LO 
QCD calculation (right), error bars on data points represent the statistical er- 
ror, the shaded area represents the systematic uncertainty. Lower: percentage 
difference between data and theory relative to theory for the dijet mass cross 
section measured in _D0 using 34 pb~^ of Run II data; the NLO theoretical cal- 
culation (JETRAD, ^ = P^"°^/2, Rsep = 1-3) involves the CTEQ6M (left) 
and MRSTOl (right) PDF sets. Data points include the statistical error, the 
solid band represents the total systematic uncertainty. 



are summarized in the upper plots of fig0| The comparison to Run I data 
(upper left) shows that the Run II dijet mass spectrum already extends the 
Run I one by about 350 GeV ^ . The effect of the increase in cross section as 
consequence of the increase in y^ is more evident in terms of the Run II over 
Run I cross section ratio (upper right). A good agreement with a LO calculation 
is observed within systematic uncertainties (dominated by the one on jet energy 



^The highest dijet mass value was Mjj = 1364 GeV. 



scale). The ongoing work performed in order to optimize the jet corrections, 
to use the MidPoint algorithm as well as to implement a comparison to a NLO 
calculation is expected to update these preliminary results. 

The lower plots of figE] summarize the results obtained from D0 rela- 
tively to an integrated luminosity of 34 pb~^. Here the dijet mass spectrum 
measurement is compared to a NLO calculation (JETRAD) on a linear scale 
considering two different choices for the PDF sets: CTEQ6 (lower left) and 
MRSTOl (lower right). As in the case of the inclusive jet cross section, a rea- 
sonable agreement between data and NLO theory is observed within the errors 
(dominated by the jet energy uncertainty) for both the PDF sets. 

4 W + Jets Production 

The production of W bosons in association with high energy hadronic jets 
at the Tevatron pp collider provides the opportunity to test pQCD at large 
momentum transfer. Understanding the QCD production for W -I- jets events is 
also important for the study of many Standard Model and new physics processes 
(such as top quark measurements and Higgs and Susy searches) for which it 
represents a major background. 

A preliminary study was performed in CDF using 127 pb^^ of Run II data 
at ^/s = 1.96 TeV collected in the period 2002-3. The production rates and 
kinematic properties of W -I- > n jets events were studied and compared to the 
predictions of an "enhanced" LO QCD Monte Carlo (MC). The well understood 
electroweak W — > et/ decay channel was considered as it provides an efhcient 
and clean identification of W candidate with a low background contamination. 
Jets, reconstructed with the standard CDF cone algorithm (JetClu) with R = 
0.4^^ and corrected to the parton level, were selected by requiring Et > 15 
GeV and |r/| < 2.4. QCD production, faking W -I- n jet events, was found to 
be the largest source of background in all jet multiplicity bins up to n = 3, 
with top production becoming dominant at higher (> 4) multiplicities. The 
inclusive cross sections for pp — > [W — *■ e v) + n jets production {(Tw>njets) 
were obtained after backgrounds and detection efficiencies (as derived from 
data and MC studies) were properly taken into account. The "enhanced" 

LO QCD MC consisted of the ALP GEN '^ program {fi = M^, <{P^f>] 
CTEQ5L PDFs), used to generate W ^ ei^ + n partons (n = 1 to 4) at LO, 

interfaced with HERWIG *^ for the implementation of the shower evolution 
for the initial and final state radiation, of the hadronization process and for 
the simulation of the undelying event; finally, the full detector simulation and 
event reconstruction were applied. 



■*The selection of this cone size was motivated by related studies on top 
quark physics. 
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Figure 5: Upper: W + > n jets cross sections measured in CDF analyzing 127 
pb^^ of Run II data at \/s ~ 1.96 TeV compared to LO theoretical predictions 
(left); ratio of cross sections (a'w>njets/(^w>n-ijets) as a function of jet multi- 
plicity (right). Error bars on data points include both statistical and systematic 
uncertainties. The filled band corresponds to the variation of the theoretical 
predictions with fiR/p- Lower: differential cross section as a function of jet Et 
for the leading jet in W > 1 jets events, the second highest jet Et in W > 2 
jets events and so on up to W > 4 jete events, respectively from top to bottom 
(left); Mjj distribution for the two leading jets in W > 2 jets events (right). 
Data points include statistical errors, while the band represents the jet energy 
systematic uncertainty. The two lines are fits to the theoretical distributions 
for two different values of ^^ip. 



The measured W + > n jets cross section as a function of the jet multipHc- 
ity is compared to the theoretical calculations in fig[Sl (upper left). Error bars 
on data points represent both statistical and systematic uncertainties which, 
being dominated by the systematic error on jet energy and ranging from ~ 



13% for aw>ijets to ~ 45% for aw>4jets, clearly limit the sensitivity of this 
measurement. The filled band shows the variation of the theoretical prediction 
with the renormalization/factorization scale {fJ-R/p): the dominant theoretical 
systematic uncertainty as expected for a LO calculation^. Also shown in figEl 
(upper right) is the ratio aw>njets / o'wyn-ijets as a function of the jet multi- 
plicity which gives a measure of the decrease in cross section with the addition 
of 1 jet. Being characterized by a reduced dependence on systematic uncer- 
tainties for both data and theory measurements this ratio is related to the 
magnitude of ag (even if not giving a direct measure of it). For comparison, 

Run I results'"'^ are also reported in the same figure. The measured Run II 
over Run I cross section ratios are found to be in agreement with theoretical 
expectations. In general, it can be concluded that a reasonable data to theory 
comparison is observed. 

Furthermore, the theory to data comparison was also investigated in some 
jet kinematic variable distributions: the differential cross section as a func- 
tion of the jet Et and the dijet invariant mass (Mjj) and angular separation 
(ARjj). Figl^l (lower left) shows the distributions for the leading jet Et in W 
> 1 jets events, the second highest jet Et in W > 2 jets events, and so on up 
to W > 4 jets events, respectively going from the top to the bottom plot. Fig|5l 
(lower right) reports the Mjj distribution for the two leading jets in W > 2 
jets events. Data points include statistical errors, while the band represents 
the jet energy systematic uncertainty (the dominant one). The two lines are 
fits to the theoretical distributions calculated at different values for ^Ir/f- In 
general, a fair data to theory comparison is found within errors^ especially for 
the Mjj and ARjj distributions (which, at low values, are in particular sen- 
sitive to soft and coUinear jet production) indicating that the LO calculation 
convoluted with the HERWIG parton shower approach, even if representing a 
partial higher order correction, can reproduce the data. However, the not neg- 
ligible discrepancies between data and theory in the differential cross section at 
high Et can indicate some limitations of the parton radiation in the shower 
approach in properly describing the high multiplicity topologies'". 



^The inclusive W cross section {aw>ojets) was generated at LO by HER- 
WIG. This LO calculation, being independent of QCD effects, is not sensitive 
to the choice for jjlr/f and is lower than data due to the LO approximation. 

® Theoretical uncertainties are characterized by a reduced dependence on 
Mfl/F with respect to the predictions on (Tw>njets- 

^A major theoretical limitation in such kind of studies comes from the dif- 
ficulty to proper merge matrix element calculations for different parton multi- 
plicities avoiding double counting. For this reason only inclusive quantities are 
considered. 



5 Underlying Event Studies 

The hard scattering process at hadron colliders is usually accompanied by 
the so called "enderlying event" (UE) which, consisting of the contributions 
from beam-beam remnants, initial and final state radiation and ("semi-hard") 
multiple parton interactions, has to be removed in order to achieve precise 
comparisons between jet measurements and pQCD predictions. Consequently, 
an accurate modeling of the UE is important for all analyses involving jets 
in the final state, especially at low energies. At the same time our current 
understanding of this precess is limited as it involves both perturbative and 
non-perturbative QCD. It is clearly important to check how current QCD MC 
models describe the observed properties of the UE and, if possible, to device 
an optimal tuning so that to improve their fitting to data results. 

Such kind of studies were performed in CDF during Run I LLLl' . The 
"transverse" region perpendicular to the leading jet of the event in the az- 
imuthal angle (see figEI left) is expected to be very sensitive to the UE and it 
was studied using charged particles ®. In general, it was observed that HER- 
WIgEIi, ISA jet US) and PYTHIaCS qcD MC models with their default 
parameters do not describe correctly all the properties of the "transverse" re- 
gion when compared to minimum bias and jet trigger data. In particular, it 
was found that a simple minimum bias modeling for the "beam-beam rem- 
nants" contribution to the UE (HERWIG, ISAJET) is not able to account for 
the observed charged particle multiplicity while, with the inclusion of multiple 
parton interactions, PYTHIA gives better results. It was actually found that 
PYTHIA, with an adequate tuning of its parameters in order to enhance this 
contribution, is able to fit data very welll^i^. 

Similar studies have been performed in CDF using Run II data at ^/s = 
1.96 TeV. As in Run I, the topological structure oipp collisions was considered 
in order to make a phenomenological study of the UE in minimum bias and 
jet trigger data samples to be compared to HERWIG and PYTHIA QCD MC 
predictions. In the study reported here, the direction of the leading calorime- 
ter jet (JetClu R = 0.7, \r]\ < 2) was used to isolate regions of r]-(j) space that 
are sensitive to the UE (figEl left) and observables related to charged particles 
(reconstructed by the central tracking system with Pt > 0.5 GeV/c and |?7| 
< 1) were considered. Two classes of events were defined according to the jet 
topology: "leading jet" events, with no restriction applied on the second high- 
est Et jet and "back-to-back" events, where the two jets were required to be 



*As the study of the UE mostly involves very low energy particles, only 
the charged particle component of the UE was considered because the CDF 
tracking system (immersed into a 1.4 T magnetic field) measures the momenta 
of low Pt tracks more accurately than the calorimeter measures their energies. 
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Figure 6: Left: The leading jet direction is used to define three regions in 
azimuthal angle (each spanning 120°^ with different hadronic activities in the 
event. The "transverse" region is very sensitive to the UE contribution. Right: 
average charged particle Pt sum density (<dSPT/drid(j)>) in the "transverse" 
region as a function of the leading jet Et (Et ) for "leading jet" and "back-to- 
back" events (as defined in the test) in Run II data, compared with predictions 
from HERWIG and PYTHIA (tuned on Run I data: "Tune A"). 



nearly back-to-back (A0i2 > 150°, Et'^'^/Et'^^ > 0.8). This classification was 
introduced in order to select a subsample ( "back-to-back" configuration) where 
hard initial and final state radiation are suppressed to increase the sensitivity 
of the "transverse" region to the "beam-beam remnants" and multiple parton 
scattering component of the UE. FigEl (right) shows the average charged par- 
ticle Pt sum density {<d'E,PT/dr]d(j)>) in the "transverse" region as a function 
of the leading jet E^ {Et'^^) in "leading jet" and "back-to-back" events. Also 
shown are the predictions from HERWIG and PYTHIA after a full detector 
simulation. The charged track activity is different for the two configurations: 
for the "leading jet" case the density rises with increasing Et , while for the 
"back-to-back" one it slightly falls with increasing Et ■ The rise for "lead- 
ing jet" events is attributed to hard initial and final state radiation which has 
been suppressed in "back-to-back" events whose opposite trend with increasing 
Et might be due to a "saturation" of the multiple parton interaction at small 
impact parameter. Such effect is expected to be included in PYTHIA (with 
multiple parton interactions included) but not in HERWIG (without multiple 
parton interactions). Indeed, in figEl (right) we see how PYTHIA tuned on Run 
I data ("Tune A" , version 6.206) fits very well the data for both configurations, 
while HERWIG looks working only at higher Et . Similar results are observed 
considering the average charged particle number density (<dN/d7yd0>) as well 
as in additional studies involving, for instance, the A4> dependence of these 



densities (dSPT/dTyd^ and dN/dyyd^) relative to the direction of the leading 
jet. 

6 Conclusions 

A very exciting and important QCD physics program is ongoing at the Tevatron 
pp collider where the increase in center-of-mass energy from 1.8 to 1.96 TeV 
and the higher statistics of Run II data samples are expected to extend Run I 
results at high Et^ . 

Some preliminary Run II results have been reported in this contribution. 
The measured inclusive jet and dijet mass cross sections are in reasonable agree- 
ment with NLO pQCD (EKS, JETRAD) within the uncertainty on jet energy 
scale in the data and on gluon PDF at high x in the theoretical predictions. The 
production of W + jets events is fairly described by an "enhanced" LO QCD 
Monte Carlo (ALPGEN + HERWIG). Data to theory comparison is Hmited by 
the systematic error on the jet energy scale and by the ^h/f dependence of LO 
calculations. The underlying event is well described by the model implemented 
in PYTHIA using parameters tuned on Run I data. 
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